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We study the Brownian motion (BM) of optically trapped graphene flakes. These orient orthogo- 
nal to the light polarization, due to the optical constants anisotropy. We explain the flake dynamics, 
measure force and torque constants and derive a full electromagnetic theory of optical trapping. The 
understanding of two dimensional BM paves the way to light-controlled manipulation and all-optical 
sorting of biological membranes and anisotropic macromolecules. 



The random motion of microscopic particles in a fluid 
was first observed in the late eighteenth century, and goes 
by the name of Brownian motion(BM)[l]. This was as- 
cribed to thermal agitation [2], leading to Einstein's pre- 
dictions of the resulting particle displacements |^]. BM is 
ubiquitous throughout physics, chemistry, biology, and 
even finance. It can be harnessed to produce directed 
motion}^. It was also suggested that thermally activated 
BM may be responsible for the movement of molecular 
motors, such as myosin and kinesin[5]. When a Brownian 
particle (BP), i.e. a particle undergoing BM, is subject 
to an external field, e.g. a confining potential, the fluid 
damps the BM and, in a high damping regime, such as 
for a BP in water, the confining potential acts as a cut-off 
to the BM dynamics. This is free for short times (high 
frequency limit), while is frozen at longer times (low fre- 
quency limit) [8]. These processes have perfect ground in 
experiments with optical traps, where a BP is held by a 
focused laser beam, i.e. an optical tweezers [9|. In this 
context, BM can be utilized to investigate the proper- 
ties of the surrounding environment [10, 11], as well as of 
the trapped particle, and for accurate calibration of the 
spring constants of the optical harmonic potential 12l.[l3|. 



Dimensionality plays a special role in nature. From 
phase transitions [M], to transport phenomena [HI, two- 
dimensional (2d) systems often exhibit a strikingly dif- 
ferent behavior [14j. Nanomaterials are an attractive 
target for optical trapping[16-18]. This can lead to 



top-down organization of composite nano- assemblies |16j . 
sub-wavelength imaging by the excitation and scan- 
ning of nano-optical probes [l^, photonic force mi- 
croscopy with increased space and force resolution [Tij. 
Graphene [HI is the prototype 2d material, and, as such, 
has unique mechanical, thermal, electronic and optical 
properties (2q|. Here, we use graphene as prototype ma- 
terial to unravel the consequences of BM in 2d. 

Graphene is dispersed by processing graphite in a 
water-surfactant solution, FigHJi. We do not use any 
functionalization nor oxidation, to retain the pristine 
electronic structure in the exfoliated monolayers [2l|423|. 
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FIG. 1: a)Scheme of water dispersion of graphene. b)TEM 
of representative flake, showing the typical honeycomb struc- 
ture, c) Number of layers per flake from TEM images, showing 
up to 60% SLG.d), Raman spectrum of an optically trapped 
flake, for 633nm trapping and excitation wavelength. 



We use di- hydroxy sodium deoxycholate surfactant. High 
resolution Transmission Electron Microscopy (HRTEM) 
shows~10-40nm flakes. FigdjD is an image of one such 
flake, with the typical honeycomb lattice. By analyz- 
ing over a hundred flakes, we find~ 60% single-layer 



(SLG), much higher than previous aqueous [22| and non- 
aqueous dispersions [21), and the remainder bi-and tri- 
layers (FigUhj.We then place 75/il in a chamber attached 
to a piezo-stage with Inm resolution. Optical trapping 
is obtained by focusing a near-infrared (NIK) (830nm) 
or a Helium- Neon (633nm) laser through a 100 x oil im- 
mersion objective (NA= 1.3) in an inverted configura- 
tion. The latter is coupled to a spectrometer through 
an edge filter. This allows us to use the same laser light 
both for optical trapping and Raman scattering, realiz- 
ing a Raman optical tweezers (ROT) to directly probe 



2 




FIG. 2: a) Experimental setup. A laser beam is expanded 
to over- fill the back aperture of a high NA lens. Geometry, 
relevant angles and axes are also shown, b) Laser is switched 
on and the flake is drawn into the optical trap, c) Laser is 
switched off and the flake released. 

the structure of the trapped flake, FigHJi. In both se- 
tups, the particles are imaged through the same objec- 
tive (Figl2^) that focuses the trapping light onto a CCD 
camera, with diffraction limited resolution. FiglSt shows 
an flake drawn into the optical trap when the laser is 
switched on. When the laser is switched off (FiglSt) the 
flake is released and diffuses from the trap region. In 
both set-ups, the minimum trap power is~ 1 — 2mW. 

A typical Raman spectrum of trapped flakes measured 
at 633nm is plotted in Fig(l]i. Besides the G and 2D 
peaks, this has significant D and D' intensities, and the 
combination mode D-hD'~2950cm~^. The G peak cor- 
responds to the E2g phonon at the Brillouin zone centre. 
The D peak is due to the breathing modes of sp^ rings 
and requires a defect for its activation by double reso- 
nance (DR)!!!-!!^. The 2D peak is the second order 
of the D peak. This is a single line in SLG, whereas it 
splits in four in bi-layer graphene, reflecting the evolu- 
tion of the band structure [2^. The 2D peak is always 
seen, even when no D peak is present, since no defects 
are required for the activation of two phonons with the 
same momentum, one backscattering from the other. DR 
can also happen intra- valley, i.e. connecting two points 
belonging to the same cone around K or K'. This gives 
rise to the D' peak. The 2D' is the second order of the 
D' peak. The large D peak in Fig(T]i is not due to a large 
amount of structural defects, otherwise it would be much 
broader, and G, D' would mergej25|- We assign it to the 
edges of our sub-micron flakes [27|. We note that the 2D 
band, although broader than pristine graphene, is still fit 
by a lorentzian. Thus, even if the flakes are multi-layers, 
they are electronically almost decoupled (isj. 

Positional and angular displacements of a flake in the 
optical trap are detected by back focal plane (BFP) in- 
terferometry using the forward scattered light from the 



trapped particle [l2[ [HI. The BFP interference pattern 
is determined by the flake orientation through the rele- 
vant angles (0, 6). Since the trapped flake is aligned with 
the yz plane, fluctuations occur in the small angle limit, 
(/) <C 1,^ <C 1, and the particle tracking signals in the 
cartesian directions are: Sx ~ /^^^ (X — a0 + b^) ; Sy ~ 
Py {Y -\- c(j)) ; Sz ~ PzZ^ where Pi are detector calibra- 
tion factors, X, F, Z center-of-mass coordinates, a,b,c de- 
pend on flake geometry and optical constants. Sz is not 
affected by angular motion for small angles. 

FigEk visualizes the BM of a graphene flake, com- 
pared with a nanotube bundle (FiglSb) and a spherical 
latex microbead (FiglSh), measured in the same exper- 
imental conditions. It is clear that these 2d, Id and 
3d objects exhibit distinct behaviors. The difference in 
the dynamics is due to the particle shape and optical 
properties. For a spherical particle the hydrodynamics 
is isotropic. Therefore, the different extent of fluctua- 
tions from equilibrium is only due to the anisotropy of 
the optical potential [HI, [29|. For a linear nanostructure 
anisotropic hydrodynamics leads to a much increased mo- 
bility along the optical axisfisj. In contrast, the graphene 
flake has increased fluctuations in both longitudinal and 
transverse directions, that we ascribe to a higher contri- 
bution from rotational motion with respect to nanotubes. 
As discussed later, this is a fingerprint of the 2d geome- 
try, yielding an increased sensitivity to angular fluctua- 
tions about the optical axis. The large optical anisotropy 
of graphene enhances this further, aligning the flake or- 
thogonal to the light polarization. The effect of the rota- 
tional BM is illustrated in the histograms of Figsl3ll,e,f, 
where the contribution from a superposition of transla- 
tional and rotational fluctuations is seen in the transverse 
directions, but is absent in the longitudinal. 

In order to extract quantitative data, we first ana- 
lyze the flakes hydrodynamics, which encompasses trans- 
lational and rotational motions. The viscous drag 
and torque are described by the anisotropic mobility 
tensors [y, 0] P^^-, for translations and P^^-, for rotations. 
These are related to the fluid dynamical viscosity r] 
(0.91 ImPa s for water at 24° C) and particle size. We 
approximate graphene flakes as extremely flat ellipsoids, 
with transverse size A much larger than their height h. 
This allows us to exploit the analytic solutions for uniax- 
ial ellipsoids [6]. The hydrodynamic mobilities are then 
only function of r] and A: Py ^ 8^,r^ ^ T^A'^"^ ^ 
^^^3 . The rotational mobility is the same for any axis 
through the center-of-mass, while the parallel transla- 
tion mobility is 2/3 the perpendicular one. From TEM, 
A - 25nm. Thus Py - 5.49/im/(fN s), P^ - 8.26/im/(fN 
s), P^ - 52.6/(fN nm s). 

We describe the Brownian dynamics of trapped 
graphene using uncoupled Langevin equations 0: 

ftX,(t) = -cj,X,(t)+^,(t), i = x,y,z (1) 
dtHt) = -^Mt)^^^{t), (2) 
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FIG. 3: a) Three-dimensional BM of a flake as compared to 
that of b) a nanotube bundle, and c) a latex microbead. 10"^ 
data-points are extracted from the Si{t). d,e) Histograms of 
the transverse signals Sx{t), Sy{t). In the transverse direc- 
tion, both translational and angular fluctuations are super- 
posed. The difference in the root mean square widths of the 
fluctuations in x and y arises from the flake shape, optical 
anisotropy, and different curvatures of the optical potential 
in the directions parallel and perpendicular to the initial po- 
larization. For each graph the QPD voltage-to-position cali- 
bration factors I3i are obtained using the calculated mobility 
coefficients and amplitude of the signals' autocorrelation func- 
tions Cii{0) — P^ksT /ki for the position fluctuation contribu- 
tions only. The root mean squares transverse displacements 
from a gaussian fit: \/~(^^ = 57 it 2nm; ^JJ^ = 53 =b 2nm. 
For rotations: (02) = 0.11 dz 0.01rad'-^7deg, consistent with 
the small angle approximation.f) Longitudinal signal Sz{t), 
only due to center of mass fiuctuations;'^/ {z"^) = 217 =1= 5nm 



(3) 



where Xi^(j)^0 are stochastic variables associated with 
position and angular coordinates, ^i(t) random noise 
sources with zero mean and variance {S,i{t)£,i{t + r)) = 
2kBTTiS{r)^ uji = F^/ci, = T^k^ and Qq = T^ko re- 
laxation frequencies related to force and torque constants 
and mobility tensor components. The confining potential 
torque on the lab axes is only relevant for orientational 
dynamics, while not affecting center-of-mass motion in a 



small angle regime. Also, due to the strong yz- alignment, 
and because angular fluctuations are small, the radiation 
torque along y (aflFecting 0) is small and ^ 0. 

We now evaluate the temporal correlations between 
the particle tracking signals, which yield the trap 
parameters [Tlj. For a non-spherical particle, correla- 
tion function analysis reveals information about center- 
of-mass and angular fluctuations, hence on trap force and 
torque constants [HI . For a strongly aligned 2d flake, the 
autocorrelation of the transverse tracking signal Cu (r) = 
{Si{t)Si{t + r)) decays with lag time r as a double ex- 
ponential corresponding to positional and angular relax- 
ation frequencies uoi {i = x,?/), , whereas, since the 
stochastic variables are uncorrelated in the small angle 
regime, the cross-correlations Cxyir) = {Sx{t)Sy{t + r)) 
of the transverse signals decay as a single exponential, 
with a relaxation rate corresponding to ^l^. 

These allow us to derive the optical force constants 
from the relaxation frequency measurements. FigH^ is 
a representative autocorrelation function analysis of the 
transverse tracking signals Cii{r) = {Si{t)Si{t -\- r)) {i = 
X, y). These data are well fltted by two exponentials with 
ujx = {8.6±0.2)xl0h-\ujy = (12.9±0.3)xlOVi forthe 
translational decay rates and = (3.0 ±0.1) x 10^s~^ 
for the angular ones (obtained as average values from 
the Cxx and Cyy slow relaxation rate). FigHJ) shows 
that the autocorrelation of the axial signal Czz is well 
fitted by a single exponential with Uz = (7.70 ± 0.05) x 
lO^s-^ In Figgt the cross-correlation of the trans- 
verse signals Cxy is fitted by a single exponential with 
= (2.90 ± 0.05) X lO^s"-^, consistent with the value 
obtained from the autocorrelation functions. Repeat- 
ing these measurements over ten diflFerent fiakes, and 
using our estimation of the hydrodynamic mobility pa- 
rameters, we obtain the spring constants ki = oJi/Ti 
to be kx = 1.1 ± 0.4pN//im, ky = 1.3 ± 0.5pN//im, 
kz = 0.08 ± 0.03pN//im and torque constant about the 
propagation direction /c^ = l^^/F^ = 9 ± 3fmN-nm/rad, 
where the uncertainty takes into account the 40% spread 
on fiake size. Note how the measured force constants only 
depend on the fiake transverse size A, not on thickness 
/i, because of the 2d geometry that strongly eflFects both 
hydrodynamics and radiation force and torque. 

We calculate the radiation force Fnad and torque FRad 
within the transition matrix approach 29t 30| . We first 
consider the incident field in the focal region of a high 
NA lens[29|. FRad and FRad are derived from the lin- 
ear and angular momentum conservation for field and 
graphene [30|. The dielectric constant of graphene is 
highly anisotropic (3ll - [33| with components e± and e\\ in 
the directions perpendicular and parallel to the c axis 
(see Figl2^)[32j. At 830nm the graphene refractive in- 
dex is = 3 + il.5; ny = 1.694[32]. Note that the 
n^ imaginary part yields a large absorption, while for 
ny this is negligible. We then calculate FRad(r), r be- 
ing the fiake center of mass relative to the focal point 
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FIG. 4: a)Transverse (x,y), b) longitudinal (z) autocorrela- 
tions. The double exponential decay in a) is due to the dif- 
ferent time scales of translational and angular BM. A single 
exponential is seen in the axial (z) direction. Solid lines are ex- 
ponential fit s.c) Transverse signal cross-correlations revealing 
only angular fluctuations, having decay rate consistent with 
the slower part of the transverse auto-correlation. The nega- 
tive sign is related to the 2d geometry, resulting in an oppo- 
site x,y phase during rotation. d,e,f) Calculated optical torque 
components. Stable orientation arises for a flake in the yz- 
plane (orthogonal to the polarization axis) . For Ty and Tz no 
data is shown for >20° because the flake is expelled from 
the trap by radiation pressure. The polarization torque {Tz) is 
two orders of magnitude higher than the other components. g) 
The optical trapping efficiency components Qi — cFRad,i/^i^ 
(c velocity of light, n=1.33 water refractive index and P laser 
power) are proportional to the optical force. For small dis- 
placements the force follows Hooke's law and the derivatives 
at the equilibrium position define kx,ky, kz 



(Fig 12^). Trapping occurs when FRad(r) vanishes with a 
negative derivative, FigH^. Due to the flake and field 
symmetry in the focal region 0, trapping occurs on 
the optical axis. For small displacements from equilib- 
rium, the trap is well approximated by an harmonic po- 
tential V{xi) = ^^i^xyz^i^h wi^^ spring constants 
kz <kx^ ky depending on both fla ke g eometry and gaus- 
sian beam power and polarization (l8[[29j . 

The flake orientation is specified by the angles (/:?, 
Figl2^. For each orientation we determine the center of 
mass trapping position, then the torque relative to each 
axis at that position (Figsll]i,e,f): T = ^^s^p^FRad, 
where is the incident field amplitude, ax the flake 
extinction cross section, n = 1.33 the water refrac- 
tive index. Orientational stability occurs when T van- 
ishes with negative derivative with respect to both i?, 
(Figsl4li,e,f). We get that stable trapping is achieved 
when the flake plane is parallel to yz. When the polar- 
ization axis lies on the flake plane (e.g. when the flake is 



parallel to xy or xz)^ the radiation pressure is so strong 
that the flake is pushed out of the trap. This is a con- 
sequence of the large imaginary part of e^. As shown in 
Figsl4li,e, the flake is stable under small angle rotations 
around its equilibrium orientation, while for larger values 
of (j) and (>20°) it is expelled from the trap by radia- 
tion pressure. Moreover, the polarization torque (T^) is a 
hundred times larger than the other components because 
of optical constants anisotropy. 

The graphene dimensionality and strong anisotropy in 
both optical and hydrodynamic properties determine the 
flake stability. The ability to discriminate linear and an- 
gular motion is key to understand optical trapping of pla- 
nar structures. Our results and methodology are generic, 
and can be extended to any 2d structure, such as bi- 
ological membranes. These contain micro-domains that 
play a role in signal transduction and trafficking[34'|. The 
diffusion and BM of planar structures in a restricted, 
anisotropic environment could allow to understand the 
dynamics and function of these domains. The optical 
trap is an ideal environment for spectroscopic and me- 
chanical probing of such structures, linking their BM dy- 
namics to their form, interactions and ultimate function. 
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